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Abstract 

We derive an upper limit on the neutrino mass from the WMAP three-year data alone by 
employing a deterministic minimisation method based on a grid search in multi-parameter space. 
Assuming the flat ACDM model with power-law adiabatic perturbations, we find ^ m^ < 2.0 eV 
in agreement with the result of the WMAP team. This result, the limit being nearly the same 
as that from the WMAP first-year data, means that the fundamental limit on the neutrino mass 
obtainable from the cosmic microwave background alone is already nearly met, as anticipated from 
the previous analysis. We also clarify the role of the polarisation data in deriving the limit on the 
neutrino mass. 
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In Pi it was shown that a vahd hmit on the neutrino mass can be derived from the cosmic 
microwave background (CMB) experiment alone. This contrasted to the then beheved view 
that such a hmit is obtained only when the CMB data are combined with those from large- 
scale clustering of galaxies the limit given by Tegmark et al [2| derived from the first 
year data of the Wilkinson Microwave Anisotropy Probe (WMAP-1) ^ is Yl'^^i^ < 
at a 95% confidence. This implies that massive neutrinos alone may constitute the entire 
dark matter of the Universe. This conclusion was supported by the analysis of Elgar0y and 
Lahav 0|. Our limit from WMAP-1 alone, ^m^< 2.0 eV (95% CL) obtained from a 
deterministic grid search, contradicts with that of \2l- We ascribed this discrepancy to the 
possibility that Markov chain Monte Carlo used in [2] was not long enough to sample the 
likelihood function in the presence of strong parameter degeneracy. 

The argument given by is based on the well-known fact that massive neutrinos dimin- 
ishes the small-scale power by free-streaming and this leads to a limit on the neutrino mass. 
In this argument the fact was not taken into account that massive neutrinos affect the CMB 
perturbations not only by the modification of the power spectrum but also in some not so 
obvious ways: characteristically, the acoustic length scale and the relative heights of peaks 
are modified. The shift of the position of the acoustic peaks and the height of the first peak 
alone can be absorbed into the change of cosmological parameters, but when the heights of 
the second and third peaks relative to the first are modified, the shifts cannot be absorbed 
into that of the cosmological parameters, leading to a limit on the neutrino mass. This 
means that a limit on neutrino mass could be obtained from the observed CMB multipoles 
only if neutrinos are nonrelativistic at the recombination epoch, i.e^ ''^v-> 0-6 eV. Several 
numerical analyses made recently seem to support this conclusion p, lalZj. 

At the data release of the WMAP team for their three year observations (WMAP-3) 
they gave a limit on the neutrino mass 'Y^m^ < 2.0 eV (at a 95% confidence) from 

their CMB data alone, using a Markov chain Monte Carlo analysis. This agrees with our 
(1 

conclusion in 1], where it was shown that the limit from WMAP-1, ^m^j < 2.0 eV, could 
hardly be improved even with the increasing quality of the CMB data. In this report we 
attempt to verify this explicitly by repeating the previous analysis with the WMAP-3 data 
using the deterministic grid search algorithm. Incidentally, we try to elucidate the role of 
the polarisation data if they could tighten the constraint on the neutrino mass. 

We take the flat ACDM model with adiabatic perturbations of the power law spectrum. 
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FIG. 1: The curve of x^iujy) as a function of the neutrino mass summed over generations. The soUd 
curve uses all available data from WMAP including polarisations. The dotted curve is a constraint 
from only temperature-temperature correlations. 

The model contains 6 cosmological parameters, the mass density parameters {ui = VLiK^) for 
cold dark matter, baryonic matter and massive neutrinos, the normalised Hubble constant /i, 
the reionization optical depth r, the scalar spectral index of primordial mass perturbations 
Us and their amplitude A. The neutrino mass density parameter ujy is related to neutrino 
masses by u^, = ^my/(94.1 eV) (the Fermi distribution is assumed) and we assume three 
generations of massive neutrinos with a degenerate mass. Theoretical CMB power spectrum 
is calculated using the CMBFAST code P and by the likelihood code of the WMAP three- 
year data release jlfll . 11, 3|. We calculate the function for a given u^, by minimizing it 
over the 6 parameters. The minimization is carried out with the grid search using the Brent 
method [l^ extended to multi-parameters, as described in [l|. 
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The result is shown in Fig. The logarithm of the likelihood —2 In L for the power 
spectra including all TT, TE, EE and BB with the pixel-based method for low multipoles is 
shown by the solid line. The Xefr 3162 degrees of freedom at the vanishing neutrino mass 
is 1.037 for the parameters = 0.0221, cu„ = 0.127, h = 0.725, r = 0.091 and ra, = 0.957 
in agreement with the WMAP-3 analysis. 

The limit corresponding to A^^ = 4 (measured from the minimum) is < 0.024 or 
^ nil, < 2.3 eV. We obtain the upper limits at a 95% confidence by an integration of the 
likelihood function C = exp{— Ax^/2}, 

cj^ < 0.0215 or ^m^<2.0eV, (1) 

which corresponds to m,j < 0.67 eV. A slight difference between the limit from Ax^ = 4 and 
likelihood integral is due to the curve flatter than quadratic at small neutrino masses. 
The actual minimum of occurs at a non-zero neutrino mass ~ 1.3 eV, but 

relative to the vanishing neutrino mass is less than one, meaning that the preference of a 
finite neutrino mass is insignificant. 

The limit we obtained agrees with that from the WMAP team, endorsing the validity 
of their Markov chain Monte Carlo analysis. We may also compare eq. (0) with that from 
WMAP-1 Q|: cUj, < 0.023. The limit changes little in spite of the significantly increased 
accuracy of the CMB data, which agrees with the analysis given in l|. This is close to 
the fundamental limit that could be obtained from the CMB alone even if the data quality 
would be higher. A further improvement of the limit needs external constraints, such as 
those from galaxy clustering. 

We now briefly discuss the role of the polarisation data in deriving the limit on the 
neutrino mass. We show in Figure 1 above the x^ curve for the TT correlation only (the 
dotted curve), which gives uo^ < 0.0250 or 'Y^niy < 2.4 eV, slightly weaker than eq. (P). 
We note, however, that there is no particular sensitivity in the TE or EE correlation to the 
inclusion of the neutrino mass beyond that of the TT correlation data. The left column 
of the six panel figure. Figure |H shows the TT, TE and EE correlators for = and 
0.02 with the other cosmological parameters fixed. In the right-column panels we rescaled 
the abscissa in the way that accounts for the modification of the acoustic length scale by a 
finite neutrino mass. The two curves almost coincide for all correlators, TT, TE and EE. 
The WMAP data for TT are accurate enough to marginally distinguish the two curves at 
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FIG. 2: TT, TE and EE correlators for three cases: cuu =0 (the sohd hne) and 0.02 (the dashed 
hne) with the other parameters fixed to the WMAP mean values. In the right panels the abscissa 
is rescaled according to the modification of the acoustic length that varies with the neutrino mass. 

approximately a 95% confidence level. The data involving polarisation, however, do not 
reach this accuracy, so that they do not give new information as to the neutrino mass. We 
do not expect any effect of the neutrino mass that modifies polarisation other than those 
through the scales in acoustic dynamics. The slight improvement on the limit of the neutrino 
mass upon the inclusion of the polarisation data, as we saw above, arises from the tightened 
constraint on r. There is a significant negative correlation between uJi, and r (see Figure 
2d in 1]). The absence of the polarisation data makes r more uncertain and rather drives 
the minimum to t ^ 0, which in turn pushes the neutrino mass towards a larger value. 
In this way, the polarisation data serve to slightly improve the limit on the neutrino mass 
through a constraint on the reionisation optical depth. 

Our final remark is that the limit from the CMB data alone is the most robust result in 
view of the uncertainties in the current large-scale galaxy clustering data. They drive the 
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best fit parameter of tfie matter density eitlier way depending upon wlietlier one takes tlie 



SDSS or 2dFGRS data 



tlie limit on tlie neutrino mass is also sensitive to the choice 



of the galaxy clustering data between the two Q] as well as how much weight is given to 
these data. To go beyond the present neutrino mass limit it would be imperative to combine 
external data, but this needs proper understanding of systematic errors involved in them. 
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